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SUMMARY 


FRANOPP, FRamework for ANalysis and Optimization Problems, is a software 
aid for the study and solution of design (optimization) problems. FRANOPP pro- 
vides the driving program and plotting capability of a user generated program- 
ming system. In addition to FRANOPP, the programming system also requires the 
optimization code CONMIN (created and maintained at the NASA/Ames Research 
Center), and two user supplied codes, one for analysis and one for output. Since 
the user supplies the analysis module which defines the design problem (more 
specifically, evaluates the objective function, constraints, and optionally 
gradients for given values of the design variables), the generated system is 
extremely flexible. 

FRANOPP provides the user with five options for studying a design problem. 
Three of the available options utilize the plotting capability supplied by 
FRANOPP and provide an indepth study of the design problem at hand. The study 
can be focused on a history of the optimization process or on the interaction 
of variables within the design problem. 

The optimization module, CONMIN, is a FORTRAN module for the minimization 
of a multi-variable function subject to a set of inequality constraints. The 
optimization problem may be linear or non-linear, constrained or unconstrained. 
The basic algorithm used for the constrained problem is the Method of Feasible 
Directions. For unconstrained problems, the Conjugate Direction Method of 
Fletcher and Reeves is used. 

This report describes FRANOPP 's design and capabilities, and how to build 
and execute a programming system which incorporates FRANOPP. It is not intended 
as a stand-alone user's guide. A familiarity with CONMIN and the Control Data 
Corporation Network Operating System (NOS) is assumed. FRANOPP was developed 
and tested on the NASA/Langley Research Center (LaRC) computers using available 
software utilities. 
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INTRODUCTION 


Various algorithms exist for the solution of design (optimization) pro- 
blems, and, correspondingly, many computer codes have arisen. In general, the 
problem is to determine the values for a vector of variables, termed design 
variables, that minimizes (or maximizes) some function of those variables, 
termed the objective function, while satisfying various linear and non-linear 
constraint equations. When using a generalized code, the user requires some means 
for defining a particular problem. One method is a user supplied code which 
contains the problem dependent formulations; typically the objective function, 
constraints and, optionally, their gradients. 

FRANOPP, FRamework for ANalysis and Optimization Problems, is a software 
aid for the study and solution of design problems. FRANOPP provides the driving 
program and plotting capability for a user generated software system which 
also includes the general purpose optimization code CONMIN (created and main- 
tained by NASA/ Ames Research Center) and two user supplied codes, one for analy- 
sis (definition of the design problem) and one for output. 

FRANOPP provides the user with five options for studying a design problem. 

The user may soley execute the analysis module, request an optimization run, 
or utilize one of the three plotting capabilities supplied by FRANOPP. The 
plotting options provide an indepth study of the problem and can be focused 
on a history of the optimization process or on the interaction of variables 
within the design problem. 

A software system incorporating FRANOPP is easy to build and use. The 
user only needs to be concerned with the problem dependent aspects of the gen- 
erated programming system. The modular design of the user-generated system 
and the code within FRANOPP make the system flexible for modifications and 
additions. FRANOPP may be used as an initial testbed for large problems before 
proceeding with a more complex optimization system. 
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CAPABILITIES AND OPTIONS 


FRANOPP has five options available to the user. These are: 

ANALIZE 
PLOT LINEAR 
PLOT CONTOUR 
OPTIMIZE 

OPTIMIZE WITH PLOTS 

A general description of the capabilities of each is given below. 

Option: ANALIZE 

This option executes the user's analysis routine for a single set of design 
variable values and outputs the objective function and constraint values. 

Option: PLOT LINEAR 

The PLOT LINEAR option allows the user to examine the dependence of a spe- 
cified output quantity on a specified input variable. In the general case, the 
user defines a design variable as the input parameter to vary over a given scale 
while all other input parameters remain unchanged. Also specified is the output 
quantity which the user wishes to study, usually the objective function or a con- 
straint. This option loops through the analysis routine varying the specified 
parameter and saving both the input values and the specified output values to be 
plotted in a linear, two-variable plot. 

Option: PLOT CONTOUR 

PLOT CONTOUR, similar to the PLOT LINEAR option, allows the user to select 
an output quantity for observation while varying any two input parameters. The 
results are displayed in a contour plot. 
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Option: OPTIMIZE 


The OPTIMIZE option uses the general purpose optimization code, CONMIN 
(ref. 1) to optimize the design problem defined in the user supplied analysis 
routine. The user begins the process by supplying to CONMIN an initial value for 
each design variable and some other required parameters. Using these initial 
values, the analysis routine evaluates the constraints and objective function. 
CONMIN is entered again and produces new values for the design variables. These 
new values are passed to the analysis routine and the model is re-evaluated. The 
process continues until either a feasible design is obtained with the objective 
function minimized or the user supplied value for the limit of iterations is 
reached. 


Option: OPTIMIZE WITH PLOTS 

This option is an extension of the OPTIMIZE option. In addition to pro- 
viding the results from the optimization process, this option saves the values 
of the objective function, design variables^ and constraints at each 
iteration, thereby providing a history of the process. The results are then 
displayed in a series of plots. There are two user defined input values which 
determine whether or not the objective function and design variable plots are 
plotted as they appear, or are normalized by the initial values before plotting. 
The user has no option for the constraint plots. These are always plotted in 
the normalized form Ggct^Sgii - 1, where 6^^^ is the actual behavior 
variable value and G,,. is the allowable behavior variable value. ^ 


^ When CONMIN computes the necessary gradients for the optimization 
process (instead of the user supplied analysis module), the values saved for the 
first design variable, except at the initial and final data points, will include 
an increment for the finite difference method. Since the increments are normally 
small, they do not change the character of the plot. 
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A behavior variable is any output from the analysis, such as stresses 
or forces. 
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PROGRAMMING DESIGN 


FRANOPP provides the framework for a programming system consisting of two 
programs. The first calculates the numerical results for the five options and 
the second handles the generation of plots. Since the plotting capability does 
not apply to all of the options, the second program is not always executed. 

The input/output files and the programs are described below. The flow of exe- 
cution through the system for all options is shown in figures 1 and 2. A more 
detailed description and a listing of each FRANOPP supplied subroutine is pro- 
vided in Appendix A. 


Input/Output Files 


The input/ output files accessed by FRANOPP 's programming system are listed 
with their major attributes. The order of their appearance reflects their posi- 
tion on the program header card. 


(1) INPUT; equated to TAPES; input to program one; input to program 

two; contains all the user defined input to execute FRANOPP; 
a detailed description of contents is presented in the sec- 
tion on execution. 


(2) TAPES: output from program one; input to program two; contains 

information regarding the type of plots requested and other 
required items. 


(3) TAPE?: output from program one; input to program two; contains the 

data arrays to be plotted. 

(4) OUTPUT: equated to TAPES; output from program one; output from pro- 

gram two; contains the results of the option executed. 


Program One 

Program one consists of four modules. It contains a driving module sup- 
plied by FRANOPP, the optimization code CONMIN, and two user supplied modules, 
one for analysis (definition of design problem) and one for documentation of 
the output. An explanation of each module follows: 
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Module FAOPS. - This module drives the entire system. It contains the 
logic for all the options and provides the interface between the modules. The 
main routine FAOP reads an input command and calls the appropriate subroutine 
to execute the requested option. The subroutines are listed below with the 
corresponding option (s) in parenthesis. 

Subroutine SUBANAL (ANALIZE) 

Subroutine PLTLIN (PLOT LINEAR) 

Subroutine PLTCON (PLOT CONTOUR) 

Subroutine OPT (OPTIMIZE and OPTIMIZE WITH PLOTS) 

The subroutine layout of this module reflects the independence and separation 
of the options. The separation allows for ease in modifying and extending 
FRANOPP's capabilities. A new option, for example, can be added without regard 
for those that already exist. 

Module CONMIN. - CONMIN is a group of FORTRAN subroutines for the minimize 
tion of a multi-variable function subject to a set of inequality constraints. 
CONMIN solves both linear and nonlinear, and constrained or unconstrained pro- 
blems. The basic algorithm used for the constrained problem is the Method of 
Feasible Directions. For unconstrained problems, the Conjugate Direction 
Method of Fletcher and Reeves is used. For information concerning CONMIN 's 
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organization and parameter definitions see reference 1. 

Module ANALYS. - ANALYS is a user supplied module that evaluates the 
objective function, constraints and, optionally, their gradients, given a set 
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CONMIN has been designed so that the user can supply the gradients 
of the objective function and constraints in the analysis module or CONMIN 
will compute the necessary gradients using finite difference. When CONMIN 
computes the gradients, the design variables will be incremented by either 
absolute or relative values which are governed by the input parameters FDCHM 
and FDCH, respectively. The larger increment prevails. It is advisable to 
adjust these parameters so that the relative increment prevails in the anti- 
cipated range of the design. 



of values for the design variables. For compatibility with CONMIN, the design 
problem must be defined as a minimization problem with the constraints com- 
puted as (actual behavior variable value/allowable behavior variable value -1). 

This module may vary slightly depending on the option executed. For 
example, when executing the PLOT LINEAR and PLOT CONTOUR options, the user may 
choose to plot variables other than the design variables, objective function 
or constraints. In this case, an additional common block must be added to the 
analysis module so FRANOPP has access to the user defined variables. The 
implementation of these variations are explained in the execution section. 

Module DEFINES. - DEFINES is a user supplied module for documentation of 
the computer printout. It allows the user to define a particular design problem 
in the output by printing appropriate explanations and definitions. This 
module should output text definitions of the objective function, design vari- 
ables, constraints and any other information necessary for understanding the 
problem. The output from this module is placed directly under the title, before 
any results are printed. 


Program Two 

The second program, entirely supplied by FRANOPP, produces the plot vector 
file^ for the plotter and therefore only needs to be executed when a plotting 
option is selected by the user. Two of the output files generated by the first 
program (TAPES and TAPE7) are used as Input to program two. A brief description 
of the main program and major subroutines is given below. 

Main Program PLOT. - The main program reads the parameter IPLOT from TAPES 
to determine which set of plots was requested and calls the necessary subrou- 
tine to do the plotting. 

Subroutine LINPLOT. - This subroutine retrieves the plotting data from 
TAPE7 and calls the necessary library routines to produce the two-variable 
linear plot specified by the PLOT LINEAR option. 


^ The plot vector file contains the plotting commands used to control 
the plotter. It is created by FORTRAN calls to the LaRC graphics library. 
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Subroutine CONPLOT. - Similar to LINPLOT, this subroutine retrieves the 
plotting data from TAPE? and produces a contour plot of any output quantity 
the user chooses, with respect to any two input variables. 

Subroutine HISPLOT. - This subroutine plots the history of the optimiza- 
tion process resulting from the OPTIMIZE WITH PLOTS option. For the objective 
function and design variables, the user can request that the plots be normalized 
by the respective initial values. 

EXECUTION 

To build and execute a system which utilizes FRANOPP, the user must perform 
the following steps: 

(1) Code the user supplied ANALYS and DEFINES modules. 

(2) Modify the dimension statements in the code supplied by 
FRANOPP to the actual problem dimensions. 

(3) Create the input file. 

C4) Create the job stream. 

Each step is described in detail below and an example is presented in Appendix B. 

Step 1: User Supplied Modules 

ANALYS. - The analysis module, ANALYS, consists of a user supplied sub- 
routine(s) which evaluates the objective function, constraints and, optionally, 
their gradients. Its design and complexity depends almost entirely on the 
problem at hand and the user. To integrate the analysis module into FRANOPP* s 
programni ng system, the user must follow the seven guidelines listed below: 

(1) The analysis module must contain a subroutine named ANALYS, which, 
if desired, can call other user supplied subroutines. 

(2) There are two required common blocks which are used to input and 
output the required data. The design variables are input to the 
module through the common block /VARABLE/. The output consists 

9 



of the objective function and constraints contained in the cormon 
block /CONSTRT/. These common blocks are defined below. 

COMMON /VARABLE/ NDV parameters representing design variables. 

COMMON /CONSTRT/ objective function, NCON parameters representing 
constraints. 

The output, stored in common blocks, passes from the ANALYS module 
to CONMIN without any intermediate changes. Therefore, the ANALYS 
module must compute the values in a form acceptable to CONMIN. 

CONMIN requires the objective function to be a minimization function. 
In addition, the constraints must be computed as 


G = 


^act^^all “ 
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where is the actual behavior variable value and Gg-ji is the 
allowable behavior variable value. In this form, the constraints are 


easily classified into three groups. 


If G > 0, the constraint is violated, 

if G = 0, the constraint is critical, 

and if G < 0, the constraint is satisfied. 

If the analysis module requires any user defined input, the user can 
supply the necessary READs in the analysis subroutine(s) and adjust 
the input file appropriately. The user should incorporate a test to 
insure the input is only read the first time the analysis module is 
entered and then stored in a common block. To facilitate this, a 
copinon block /ANDATA/ which contains one variable, LOOPCNT, supplies 
the current number of calls to ANALYS. For example, on the first call 
to ANALYS, LOOPCNT is equal to one; the second time the analysis 
module is entered, LOOPCNT equals two. The logic is outlined below: 


SUBROUTINE ANALYS 


. required common blocks 
COMMON /ANDATA/ LOOPCNT 

COMMON /any-name/ all variables appearing in the 
read statement(s) 


IF (LOOPCNT. NE.l) GO TO 101 
. read statement(s) 

101 CONTINUE 


RETURN 

END 

The read statement(s) should accept data from TAPES. The location 
of this data within the input file will be described in the input 
file section. 

(5) If additional output to that given by any option is desired, the 
user can supply write statements using TAPES. Since the analysis 
module is executed several times for each iteration, the user may 
want to incorporate a test so that the output is written once for 
each iteration, or only for the final results. The variable ITER 
in the common block /CMNNl/ (listed in guideline 7) contains the 
current iteration number. To insure that the output is written 
once for each iteration, the user will need to save the value 
of ITER in a separate common block and compare the saved value 
with ITER until they are unequal. For the final results, the 
variable LOOPCNT (see guideline 4) has the value -999 on the last 
call to the analysis module. 
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(6) If the user wishes to execute either the PLOT LINEAR or PLOT CONTOUR 
options and specify plotting variables other than the design vari- 
ables, objective function or constraints, an additional common block, 
/VINOUT/, must appear in the analysis module. 

COMMON /VINOUT/ vin(2),vout 

where vin(l) - if the user is supplying the input variable for 

the x-axis (INX < 0), the variable name should 
replace vin(l); otherwise, this is a dummy variable. 

vin(2) - only applies to the PLOT CONTOUR option, it is a 
dummy variable for the PLOT LINEAR option. If the 
user is supplying the input variable for the y-axis 
(I NY < 0), the variable name should replace vin(2); 
otherwise, this is a dummy variable. 

vout - if the user is supplying the output variable 

(INY <0 for PLOT LINEAR option, INZ < 0 for PLOT 
CONTOUR option) which corresponds to the y-axis or 
contours depending on the option, the variable name 
should replace vout; otherwise, vout is a dummy 
variable. 

(7) When executing the OPTIMIZE or OPTIMIZE WITH PLOTS option, the ana- 
lysis module may need to access or supply information other than that 
stored in the required common blocks /VARABLE/ and /CONSTRT/. For 
instance, a user may want to compute the required gradients in the 
analysis module instead of having CONMIN compute them. Two addi- 
tional common blocks which facilitate the data transfer are shown 
below. The variables are defined in the CONMIN manual and the dimen- 
sions are defined in STEP2, Problem Dependent Dimension Statements. 

Common /CNMNl/ DELFUN,DABFUN,FDCH,FDCHM,CT,CTMIN,CTL, 
CTLMIN,ALPHAX,AB0BJ1, THETA, OBJ ,NDV,NCON,NSIDE,IPRINT, 
NFDG,NSCAL,LINOBJ,ITMAX,ITRM,ICNDIR,IGOTO,NAC,INFO, 

INFOG, ITER 

Common /ACON/ VLB(N1 ) ,VUB(N1 ) ,SCAL(N1 ) ,DF(N1 ) ,A(N1 ,N3) , 
ISC(M2),IC(N3) 
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DEFINES. - The DEFINES module consists of a user supplied subroutine called 
DEFINE. The sole purpose of this module is to document the problem in the output 
listing. The output produced should define the design problem by including def- 
initions of the objective function, design variables, and constraints. The 
output must be written to TAPES. 

Step 2: Problem Dependent Dimension Statements 

Both FRANOPP supplied source files, FAOPS and PLOTS, contain subroutines 
with problem dependent dimension statements. The dimension statements have been 
coded with undefined variables for the array sizes which must be replaced by 
the actual values before compilatiK)n of the system. The following table gives 
the location of these dimension statements within FRANOPP. 


FILE 

ROUTINE 

VARIABLES USED IN DIMENSIONS 

FAOPS 

SUBANAL 

Nl, N2 


PLTLIN 

Nl, N2, ITMAX 


PLTCON 

Nl, N2 


OPT 

Nl, N2, N3, N4, N5 

PLOTS 

L INPLOT 

ITMAX4 


HISPLOT 

NDV, NCON, Nl, N2, ITMAX2, ITmX4 


The variables are defined below. It is important to note that when the 
OPTIMIZE or OPTIMIZE WITH PLOTS option is executed, variables N1 through N5, 
and ITMAX correspond to the same input values used for CONMIN. For the PLOT 
LINEAR option, ITMAX, the maximum number of iterations used in the dimension 
statements, must be greater than or equal to NX, which is the number of itera- 
tions input to FRANOPP. 
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VARIABLE 


DEFINITION 

NDV 


number of design variables 

NCON 

= 

number of constraints 

N1 

= 

NDV + 2 

N2 

= 

NCON + 2 * NDV 

N3 

= 

1 + user's best estimate of the maximum number of 



constraints (including side constraints) which 
will be active at any given time in the minimize' 
tion process 

N4 

= 

maximum of NDV and N3 

N5 

= 

2 * N4 

ITMAX 

= 

maximum number of iterations the user allows 

ITMAX2 

= 

ITMAX + 2 

ITMAX4 

= 

ITMAX + 4 


To modify the variables in the dimension statements for each new problem, 
the user can either edit the source code directly, or execute a preprocessor 
which does the equivalent. The preprocessor route is favorable in that only 
the input to the preprocessor needs to be modified for each new problem instead 
of the source code. Since the editing option is self explanatory, the discus- 
sion will continue with the execution of the preprocessor available at NASA/ 

LaRC called PRE.^ 

There are four files associated with PRE. These are INPUT, OUTPUT, INFILE, 
and PREFIL. The user must supply the preprocessor control cards containing the 
variables to be changed and their corresponding values on the input file, and 
the source code on the file INFILE. The program PRE generates a new source 
file called PREFIL, which contains constant dimensions and a group of assign- 
ment statements that set the variables used in the dimensions to the constants 
supplied. The listable output from the preprocessor verifying the success of 
the run appears on the output file. 


5 

Documentation for NASA/LaRC implementation of PRE is available 
from the LaRC User Support Office. 
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The preprocessor control cards are of the form: 

routine name. (variable 1 = constant, variable 2 = constant, . . ., 
variable n = constant). 

The job stream for execution of PRE along with an example of the input is 
shown in the section on the job control stream. 

Step 3: Input File 

The input file for FRANOPP consists of the input for both programs, FAOPS 
and PLOTS, separated by an end-of-record mark. The file contains a command 
specifying which option the user wishes to execute and the necessary input for 
that option. A table for each command showing the required input follows. 

OPTION: ANALIZE-INPUT 


Card Input & Format 

Description & Comments 

Routine 

ANALIZE 

command to execute option ANALIZE 

FAOP 

<title> 

Format(4A10) 

title used in output 

FAOP 

$INPUT DV=_,_,...,J 
Namelist format 

design variable values 

SUBANAL 

<user specified> 
(optional ) 

if the user has input into analysis 
routine from TAPE5, the data is 
placed here 

ANALYS 
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OPTION: PLOT LINEAR- INPUT 


Card Input & Format 

Description & Comments 

Routine 

PLOT LINEAR 

command to execute option 
PLOT LINEAR 

FAOP 

<title> 

Format(4A10) 

title used in output 

FAOP 

$INPUT DV=_,_. 

VIN=_,INX=_,XINC=_, 

NX=_,INY_$ 

Namelist format 

variables described on next 
page 

PLTLIN 

<user specified> 
(optional ) 

if the user has input into analysis 
routine from TAPES, the data is 
placed here 

ANALYS 

<end-of-record> 

denotes end of input to program one; 
beginning of input to program two 


<labelx,labely> 

Format(A10,5X,A10) 

labels for the x and y axes of 
the generated plot 

LINPLOT 
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OPTION: PLOT LINEAR- VARIABLES 


VARIABLE 

DESCRIPTION 

DV 

array containing design variable values 

VIN(l) 

if INX < 0, initial value for the user defined input; 
if INX > 0, unncessary 

INX 

> 0, index into DV array to specify which design 
variable to use for plotting; 

< 0, user defines the input variable for plotting 
in the common block /VINOUT/ in the analysis 
routine; 

= 0, error condition 

XINC 

the increment value for the input variable to be used 
in plotting 

NX 

the number of times to increment the input variable. 

I NY 

> 0, index into array to specify which constraint 
value to use for the y-axis of the plot; 

= 0, the objective function is used as the output 
variable in the plot; 

< 0, user defines the output variable for plotting 

in the common block /VINOUT/ in the user supplied 
analysis routine 
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OPTION: PLOT CONTOUR- INPUT 


Card Input & Format 

Description & Comments 

Routine 

PLOT CONTOUR 

command to execute option 
PLOT CONTOUR 

FAOP 

<tit1e> 

Fonnat(4A10) 

title used in input 

FAOP 

$INPUT DV=_,_,..._, 
VIN=_,_,INX=_,XINC=_, 
NX=_,IiiY=_,YINC=_,NY_, 
INZ=_,NCONINT=_$ 
Namelist format 

variables described on next 
page 

PLTCON 

<user specified> 
(optional ) 

if the user has input into analysis 
routine from TAPES, the data is 
placed here 

ANALYS 

<end-of-record> 

denotes end of input to program one; 
beginning of input to program two 


<labelx,labely,labelz> 

Format(A10,5X,A10,5X,A10) 

labels for the x-axis, y-axis 
and contours of the generated plot 

CONPLOT 
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OPTION: PLOT CONTOUR-VARIABLES 


VARIABLE DESCRIPTION 

DV array containing design variable values 

VIN(l) if INX < 0, initial value for the user defined 

input variable for the x-axis; 
if INX > 0, unnecessary; 

VIN(2) if INY < 0, initial value for the user defined 

input variable for y-axis; 
if INY > 0, unnecessary 

INX > 0, index into DV array to specify which design 

variable to use for the x-axis of the plot; 

< 0, user defines the input variable for the x-axis 

in the common block /VINOUT/ in the analysis 
routine; 

= 0, error condition 

XINC the increment value for the input variable for the 

x-axis 

NX the number of times to increment the input variable 

for the x-axis (must be ± 40) 

INY > 0, index into DV array to specify which design 

variable to use for the y-axis of plot; 

< 0, user defines the input variable for the y-axis 

in the common block /VINOUT/ in the analysis 
routine; 

= 0, error condition 

YINC the increment value for the y-axis input variable 

NY the number of times to increment the input vari- 

able for the y-axis (must be £ 40) 





OPTION: PLOT CONTOUR-VARIABLES (cont'd) 


VARIABLE 

DESCRIPTION 

INZ 

> 0, index into array G to specify which constraint 
to use for the contours; 

= 0, the objective function is used as the output 
variable in the plot; 

< 0, user defines the output variable for plotting 
in common block /VINOUT/ in the user supplied 
analysis routine 

NCONINT 

the number of contour intervals (default = 20) 
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OPTION: OPTIMIZE or OPTIMIZE WITH PLOTS-INPUT 


Card Input & Format 

Description & Comments 

Routine 

OPTIMIZE or 
OPTIMIZE WITH PLOTS 

commands to execute options 
OPTIMIZE and OPTIMIZE WITH PLOTS 

FAOP 

< title > 
Format(4A10) 

title used in output 

FAOP 

SCONPAR <CONMIN 
parameters> $ 

for description of CONMIN 
parameters, see ref. 1 

OPT 

<user specified > 
(optional ) 

if the user has input into analysis 
routine from TAPES, the data is 
placed here 

ANALYS 

<end-of- record > 

denotes end of input to program one; 
beginning of input to program two 


<normobj ,normdv > 
Format(2I2) 

input only for OPTIMIZE WITH PLOTS 
option; keys to determine whether 
or not to normalize the objective 
function and design variable plots; 
described on next page 

HIS PLOT 
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OPTION: OPTIMIZE or OPTIMIZE WITH PLOTS-VARIABLES 


VARIABLE DESCRIPTION 


NORMOBJ = 1, normalize the objective function values by 

the initial value before plotting; 

= 0, plot the objective function values without 
normal izati on 

NOTE: 1 - Any other input defaults to 0. 

2 - If the initial value of the objective 

function is 0, NORMOBJ will be disregarded 
and the plots will not be normalized. 

NORMDV = 1, normalize each design variable array 

by the starting value before plotting; 

= 0, plot the design variable arrays without 
normalization 

NOTE: 1 - Any other input defaults to 0. 

2 - If the starting value of any design 
variable is 0, NORMDV will be dis- 
regarded and all design variables 
will be plotted without normalization. 
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step 4: Job Control Stream 

The CDC job control stream for the execution of FRANOPP is available in a 
procedure file, FAOPPRC. This procedure file uses the preprocessor PRE to 
change the variable dimensions to constants in the FRANOPP supplied programs. 
FAOPPRC is designed to check the value of register one to determine if the 
plotting program needs to be executed. If plotting is requested (options PLOT 
LINEAR, PLOT CONTOUR, or OPTIMIZE WITH PLOTS) register one should be assigned 
the integer value one and the preprocessor input for both programs must be 
supplied in the control stream. Otherwise, register one is set to zero, and 
only the preprocessor cards for program FAOPS should be supplied. 

The procedure file has four parameters which the user must define in the 
call statement. These are: 

A1 = binary for DEFINES module 

B1 = binary for analysis module 

Cl = input file for FRANOPP 

UNI = user number where A1 , B1 , and Cl reside 

An example of a user's job control stream is given below, followed by a 
listing of the procedure file FAOPPRC. 

User's job stream 


<D£L1VERV INFOt1»»TIOM> 


/JOB 

FR ANOP . T 1 «• . Cn 1 1 0009 . 

USER, (USER N0.>. 

CHARGE , (CHARGE NO.>.LRC. 

DELIVER. (DELIVERY INF0MAT10N> 

RFL( 110000) 

SETiRiM) ^ Set register 1 to 1 for plotting 

GET. FAOPPRC/UN* (USER N0.>. 

CALL,FA0PPRC(A1-EXDEFB,B1-EXANAI.C1*EXINPT,UN1-(USER N0.>) 

Plotting instructions 

EXIT. 

/EOR 

SUIAMAL(m«6,Na-ll) 

PLTLIM(H1 

PLTCONCMl 
0PT(M1 
/EOR 

HISPl0T(Nl-S^ME*ll,NDU*H,MC0N-3,ITI1AXa*4a,ITI1AX4M4) 


iL(m«6,Na-ii) ^ 

l(Hl-G.Na-ll.ITnAX-40) \ InoiJt 

i(Mi*6,Na*in I 

-6,Na'll,N3*4,N4H,N5*8)>' 


to PRE for FAOPS 


^ Input to PRE for PLOTS 
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FAOPPRC - procedure file 


COnrCNT.ttttt lEGIN PROC FAOPPRC 
CO««£NT. TO CALLS CALL. FAOPPRC (A1 • — ,» 1 ■ 
CO«i1ENT. WHERE s 

COnflENT. A1 • DEFINE NODULE BINARV FILE 
CONNEMT. B1 • ANALYSIS NODULE E INARY FILE 
CONNENT. Cl • INPUT FILE TO FRANOP 
CORNENT. UNI • ACCOUNT NUNBER L'hERE 
CONNENT. Al.ll. AND Cl ARE 

COWIENT. LOCATED. 

RAP. OFF. 

GET,PRE/UN*738625C. 

GET.FA0PS/UN-743575C . 

PGE< . .FAOPS.PREFAOP ) 

REUIND.PREFAOP. 

FTN(I-PREFAOP,B»FAOPI,L*0) 

RETURN. FAOPS.PREFAOP. 

IFfRl'0)GOTO,l. 

GET.PL0TS/UN-743579C. 

FRE(.. PLOTS. PREPLOT) 

REWIND. PREPLOT. 

FTN<I-PREPLOT.B-PLOTI,L-0) 

RETURN, PLOTS, PREPLOT. 

1,GET,A1,B1,C1/UN-UN1. 

GET, BNEUC/UN-012G8GN . 
ATTACH.FTNnLIB/UN-LIlRARV. 
LDSET.LIB-FTNHLIB, PRESET-ZERO. 
LOADIFAOPB.Al.Bl.BNEUC) 

EXECUTE!, Cl. PLTINPT.PLTDATA) 

RETURN. A l.Bl.BNEUC. 

IF<Rl-0)GOTO.a. 

REWIND, PLT INPT , PLTD ATA . 

GET. PLOT3D/UN- 159 1 37N . 

ATTACH ( LRCGOSF/UN -LIBRARY) 

LDSET( LIB-PL0T3D/LRCG0SF , PRESET- INDEF ) 

PLOTB. C 1. PLT INPT. PLTDATA. 

a,COnHENT.»«»»» END PROC FAOPPRC ttttt 


— .C1----.UN1----) 


Execution of PRE for FAOPS 


-> Execution of PRE for PLOYS 


^Execution of FAOPS 


^ Execution of PLOiS 
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CONCLUDING REMARKS 


FRANOPP, FRamework for ANalysis and Optimization Problems, is a software 
aid for the study and solution of design (optimization) problems. It provides 
the framework for a user generated programming system in which the user needs 
to be concerned with the problem dependent aspects of the system. The modular 
design, at both the system and code levels, makes the system flexible for modi- 
fications and additions. FRANOPP may be used as an initial testbed for large 
problems before preceding with a more complex optimization system. 

FRANOPP was developed and tested on the NASA/LaRC computers using the 
available software packages. Usage at another installation may require signi- 
ficant modifications. In addition, the optimization code incorporated into 
FRANOPP uses the Method of Feasible Directions for constrained problems and 
the Conjugate Direction Method of Fletchers and Reeves for unconstrained pro- 
blems. If other optimization techniques are desired, revisions to FRANOPP will 
be necessary. 



APPENDIX A; SOURCE CODE AND DESCRIPTIONS 

This appendix contains the source code for the FRANOPP supplied programs. 
Each routine is presented separately with a descriptive outline preceding it. 
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ROUTINE NAME: FAOP 

FILE NAME; FAOPS (source) , FAOPB(binary) 

PROGRAM 1 

PROBLEM INDEPENDENT 
OPTION(S); all 

INPUT: INPUT(TAPEB) - user supplied 

CARD(1 )-CARD(5) - command to specify the option to be executed 
TITLE(1)-TITLE(4) - title for output 

OUTPUT: OUTPUT (TAPE6) 

TITLE(1)-TITLE(4) 

if an unrecognizable command is read, an error message is generated 
COMMON: /PLT/ IPL0T,TITLE(4) 

SUBROUTINE CALLS(S): SUBANAL,PLTLIN,PLTCON,OPT 

DESCRIPTION: FAOP is the main routine in program one. This routine reads 

an input command and calls the appropriate subroutine to 
execute the option requested. 
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PROGRAH FAOP< INPUT, TAPE3, TAPE?. OUTPUT, TAPE5>INPUT,TAPE6-0UTPUT) 
DINENSION CARD(5) 

COnnON /PLT/ IPL0T,TITLE<4) 

READ INPUT COnriAND AND PROBLEM TITLE 


REA0(5,9ei) (CARD(I),1«1,S) 

Ml FORMAT (5A4) 

REA0(5,i010) TITLE 
FORMAT<4A10) 

URITE(6,102e) TITLE 

1620 F0RMAT(1HI,6H »»« ,4A10,5H «»*,//) 

DEFINE DESIGN VARIABLES AND CONSTRAINTS IN OUTPUT 
CALL DEFINE 

DETERMINE WHICH OPTION IS REQUESTED 
and call THE APPROPRIATE SUBROUTINE 

IF (CARD(1).NE.4H0PTI) GO TO ICl 
IPLOT«0 

IF (CARD(4).EQ.4HH PL) lPLOT-3 
CALL OPT 
GO TO loee 

IF (CARD(1).NE.4HANAL) GO TO 200 
CALL SUBANAL 
GO TO 1000 

IF (CARD(1).NE.4HPL0T) GO TO 801 
IF (CARD(2).E0.4H LIN) CALL PLTLIN 
IF (CARD(2).EQ.4H CON) CALL PLTCON 
GO TO 1000 

FORMATUHllsTHE^FOLLO^ IS NOT VALID INPUTl»,/,5X, 

2 5A4) 

1000 STOP 
END 


101 


200 


801 

SSI 
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ROUTINE NAME: SUBANAL 

FILE NAME: FAOPS( source ) , FAOPB(binary) 

PROGRAM 1 

PROBLEM DEPENDENT; dimension statements 

OPTION(S): ANALIZE 

INPUT; INPUT(TAPE5) - user supplied 

Namelist /INPUT/ DV - array of design variable values 

OUTPUT; 0UTPUT(TAPE6) - the results of the analysis 

COMMON: /ANDATA/ LOOPCNT 
/CONSTRT/ 0BJ,G 
/VARABLE/ DV 

SUBROUTINE CALL(S): ANALYS 

DESCRIPTION; This routine contains the necessary logic for the ANALIZE 
option. It reads in the values of the design variables, 
performs an analysis, and outputs the results. 
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SUIROUTINE SUBANAt 
COnnON /UARABLE/ DU(N1 ) 

COnnON /CONSTRT/ 0BJ,G(N2) 
COnnOM /ANDATA/ LOOPCNT 
NAHELIST /INPUT/ DU 

READ UALUES OF DESIGN UARABLES 

READ(S» INPUT) 

ANALYSIS 

LOOPCNT • 1 
CALL ANALVS 

WRITE RESULTS 


901 

962 

903 

904 


URITE(6,901> 

FORNATaHl. (DESIGN UARABLESt) 
URITE(6,S02) DU 
F0RnAT(4X,E12.4) 

URITE<S,903) 

FORnAT(lH0,*OBJECTIUE FUNCTIONI) 

URITE(6,9e2) OBJ 

URITE<6,964) 

F0RNAT( 1H6, (CONSTRAINTS* ) 

WRITE (6, 962) G 

RETURN 

END 
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ROUTINE NAME: PLTLIN 

FILE NAME: FAOPS( source) , FAOPB(binary) 

PROGRAM 1 

PPROBLEM DEPENDENT: dimension statements 

OPTION(S): PLOT LINEAR 

INPUT: INPUT(TAPE5) - user supplied 

Namelist /INPUT/ DV,VIN(1 ),INX,XINC,NX,INY - includes values for 
the design variables, which two variables (one 
input and one output) are to be plotted, and 
parameters that determine the scale of the input 
variable to be plotted. 

OUTPUT: OUTPUT (TAPE6) - a listing of the values to be plotted 

TAPES and TAPE? - all the necessary data for the plotting 

COMMON: /ANDATA/ LOOPCNT 

/CONSTRT/ 0BJ,G 
/PLT/ IPL0T,TITLE(4) 

/VARABLE/ DV 
/VINOUT/ VIN(2),V0UT 

SUBROUTINE CALL(S): ANALYS 

DESCRIPTION: This routine contains the necessary logic for the PLOT LINEAR 
option. The input consists of values for the design variables, 
which two variables (one input and one output) are to be plotted 
in a linear plot, and parameters that determine the scale of the 
input variable to be plotted. The input variable is incremented 
by a user defined amount in a loop containing a call to the ana- 
lysis routine. All necessary plotting information is written to 
the data files used by the second program. 
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c 

c 

c 


M 6 

987 


C 

C 

C 

c 

c 


9d9 

tea 

103 

ICfi 


993 

101 


c 

c 

c 


992 

ICS 


SUIROUTINE PLTLIN 

DIHENSION XARRAY( ITrWX ), YARRAV( ITflAX ) 

COnnON /UARABLE/ DU(Nl) 

COnnON /CONSTRT/ 0BJ.GCN2) 

COfirtON /UINOUT/ UIN(a).UOUT 
COnnON /ANDATA/ LOOPCNT 
COnaON /PLT/ lPL0T,TITLE(-4) 

NftrSELlST /INPUT/ DU.INX.UIN.XINC.NX.IHY.UOUT 
NAHELIST /PLOT/ NX.XARRAY, VARRAV 
NAnSLIST /lOPT/ IPLOT 

READ DATA 

READ(5, INPUT) 

URITE(6.906) TITLE 
FORnAT(lHl,4A10) 

URITE(6.907) 

FOROAT ( /, 5X , BXARRAY* , 1 0X, » VARRAV*. / ) 

PERFORM ANALYSIS NX TIMES 
STORE SPECIFIED RESULTS FOR PLOTTING 
IN ARRAYS XARRAY AND VARRAV 

DO 101 I-l.NX 
LOOPCNT- I 
CALL ANALYS 

IF (IMX.GT.0) CO TO 103 
IF (INX.LT.0) CO TO 102 
l'9ITE(6,909) 

FORnAT<lHi,»lNUALID INDEX - INX • 0*,/» 

, UNO PLOTTING ATTEMPTED*) 

C-0 TO 105 
XARRAV(I)-UIN(1) 

UlfUD-UINlD+XINC 
GO TO 10S 
XARRAY(I)-OU(IHX) 

DU(INX)-DU(INX)+XINC 

VARRAVa)-OBJ 

IF aNY.GT.0) YARRAV(I)-G(INY) 

IF (INY.LT.0) YARRAY(I)-UOUT 
l!SITE(6,908) XARRAYCD.YARRAVCI ) 

FORMATaH ,E12.4,4X,E12.4) 

CONTINUE 

L'RITE PLOTTING FILES FOR PROGRAM TUO 

IPLOT- 1 
URIYEt3,I0PT) 

L'RITEl 3,932) TITLE 

FORMATKAIO) 

k!RITE(7,PL0T) 

CONTINUE 

cetu:;m 

END 
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ROUTINE NAME: PLTCON 

FILE NAME: FAOPS(souce) , FAOPB(binary) 

PROGRAM 1 

PROBLEM DEPENDENT: dimension statements 

OPTION(S): PLOT CONTOUR 

INPUT: INPUT(TAPE5) - user supplied 

Namelist /INPUT/ DV,VIN,INX,XINC,NX,INY,YINC,NY,INZ,NCONINT - 

includes values for the design variables, which 
three variables (two input and one output) are 
to be plotted, and parameters determining the 
scales of the plot 

OUTPUT: OUTPUT (TAPES) - a listing of the values to be plotted 

TAPE3 and TAPE? - all the necessary data for the plotting 

COMMON: /ANDATA/ LOOPCNT 

/CONSTRT/ 0BJ,G 
/PLT/ IPL0T,TITLE(4) 

/VARABLE/ DV 
/VINOUT/ VIN(2),V0UT 

SUBROUTINE CALL(S): ANALYS 

DESCRIPTION: This subroutine contains the code for the PLOT CONTOUR option. 

The logic is the same as that for PLTLIN except that two input 
variables are defined by the user and varied in the loop. The 
output files store the necessary data for the second program, 
which produces a contour plot. 
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c 

c 

c 


907 

909 
106 

910 

107 

C 

C 

C 


C 

C 

C 


C 

C 

C 


9t8 


101 


108 


C 

c 

c 


802 

Its 


SUBROUTINE PLTCON 
DinENSION Zl^O.^O) 

CORDON /UARABLE/ DUCNl) 

CORRON /CONSTRT/ 0BJ,G(N8) 

COnnON /UIHOUT/ UINCS), uout 
CORRON /ANDATA/ LOOPCNT 

nmeust**/input/°duIuin^inx,xinc,nx,iny,yinc,nv,uout,inz,nconint 
N jyiELIST /PLOT/ XNIN.XnAX,XINC,VniN,Vf1AX,VINC,2.NC0NINT 
NARELIST /lOPT/ IPLOT 

READ DATA 


NCONINT-80 
UINCl )*UIN(8)»0. 

U0UT*0. 

READ (5, INPUT) 

URITEl 6,907) 

FORRAT < /, 5X, tXARRAV* , lOX , *YARRAV* , lOX , »ZARRAV» , / ) 

IF (INX.NE.0) GO TO 106 

URITE(6,909) 

FORFIATdHl.JlNUALID INDEX - INX • 0t,/, 

, t NO PLOTTING ATTERPTEDt) 

GO TO 105 
XtllN-UIN(l) 

IF (INX.GT.0) XWIN-DUdNX) 

IF (INY.NE.O) GO TO 107 
L'RITE(6,910) 

F0RRAT(1HI,»INUALID INDEX - INV ■ 0»,/, 

. S NO PLOTTING ATTEMPTED*) 

GO TO 105 
vniN-uiNia) 

IF (INV.GT.0) VniN-DU(INV) 


LOOP FOR X-AXIS 


DO 102 I*1.NX 
LOOPCNT -1 

IF (INX.GT.0) UIN(1)*DU(INX) 
IF (I.EQ.NX) XriAX*UIN(l) 
UIN(2)'VniN 

IF (INY.GT.0) DU(INV)-vniN 
LOOP FOR V-AXIS 


DO 101 J-l.HY 

IF (INY.GT.O) UIN(2)-Dy(INY) 
IF (J.EQ.NV) YnAX-UIN(2) 

PERFOSn ANALYSIS NXtNY TINES 


CALL ANALYS 
2(I,J)*0BJ 

IF (INZ.GT.O) Z(I,J)-G(INZ) 

IF (IMZ.LT.O) Ztl,J)*UOUT 
URITE(6,898) U1N(1 ),U1N(2),Z(1.J) 
FC.-:RAT(1H ,Eia.4,3X.E12.4,3X,E18.4) 

UIN(a)-OIN(a)+YINC 
IF (INY.GT.O) DU(INY)*UIN(a) 

CONTINUE 

UIN(1)«UIN(1)+XINC 
IF (INX.GT.0) D0(INX)-UIN(1) 

CONTINUE 

WRITE PLOTTING FILES FOR PROGRAN TWO 


IPLOT-a 
yRITE(3,I0PT) 
WRITE! 3,902) TITLE 
FCnRAT<4A10) 

WRITE (7,PL0T) 

CONTINUE 

RETURN 

END 
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ROUTINE NAME: OPT 

FILE NAME: FAOPS ( source ) , FAOPB(binary) 

PROGRAM 1 

PROBLEM DEPENDENT: dimension statements 

OPTION(S): OPTIMIZE and OPTIMIZE WITH PLOTS 

INPUT: INPUT(TAPE5) - user supplied 

Namelist /CONPAR/ CONMIN parameters 

OUTPUT: 0UTPUT(TAPE6) - CONMIN output 

TAPES and TAPE7 - if plots requested, these files contain all the 
necessary data for the plotting 

COMMON: /ACON/ VLB(N1 ) ,VUB(N1 ) ,SCAL(N1 ) ,DF(N1 ) ,A(N1 ,N3) ,ISC(N2) ,IC(N3) 
/ANDATA/ LOOPCNT 
/CNMNl/ CONMIN parameters 
/CONSTRT/ OBJ,G 
/PLT/ IPL0T,TITLE(4) 

/VARABLE/ DV 

SUBROUTINE CALL(S): CONMIN, ANALYS 

DESCRIPTION: This subroutine contains the logic for both optimization options 
(with and without plots). It reads in the CONMIN parameters 
and makes the necessary data connections with CONMIN. The opti- 
mization is performed in an optimization loop which contains a 
call to both CONMIN and the analysis routine. If plotting is 
requested, it writes the necessary data files for the second 
program. 
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SUB90UTINE OPT 

OirSNSrON S(N1 ).Gl(N2),Ca(Na),B(N3,N3),C(N4),nSl(N5) 

DinENSION SAUXtNl ).SAUG(Na> 

IMTEGSn SAUITER 

COnnON /AGON/ ULB(N1 ),UUB(N1 ),SCAUN1),DF(N1 ),A(N1,N3). 

. ISC(N2),IC(N3) 

COnnON /CKHNl/ DELFUN,DABrUN,FDCH,FDCHn,CT,CTniN,CTL,CTLniN, 
ALPHAX, ABOB J 1 . THETA, OB J , NDU. NCON, NS 1 BE , IPR INT , 

. NFDG.NSCAL.LINOBJ, ITMAX, ITRfI, ICNDIR, IGOTO.NAC, 

. INFO, INFOG, ITER 

COnnON /OARABLE/ X(N1) 

CO~aON /CONSTRT/ AOBJ.GCNB) 

CCraON /ONDATA/ LOOPCNT 
C0r,n0N /PLT/ IPL0T,TITLE(4) 

NOnELIST /COHPAR/ INFOG, INFO, NFDG, IPRINT,NDU, ITHAX,NCON,NSIDE, 

. ICNDIR,NSCAL,FDCH,FDCHn,CT,CTniN,CTLniN, THETA, 

phi,delfun,dadfun,linobj,itr;i,x,i;lb,uub, 

. Nl,Na,N3,N4,N5,nLRHAX,AB0BJl,CTL,ISC,SCAL 

NAHELIST /PLOT/ SAUITER, SAyOBJ,SftU:<,£AUG 
NOflELIST /PLTINPT/ NOITER, ITnAX,NDU,NCCN 
NAHELIST /lOPT/ IPLOT 

READ THE PARAnETERS FOR CONfUN 

READ( 5 ,C 0 NPAR) 

t'RITE<6,C0NPAR) 

MLin=nnAxa(NDU+5) 

NON-ITERATIUE PART OF ANALYSIS 

I GOTO " 0 
SAUITER— 1 

ITERATIUE PART OF ANALYSIS 


DO 1000 I 
LCOPCHT’I 


l.NLin 


CALL THE OPTiniZATION ROUTINE C0W1IN 

CALL CONPIIN(X,ULB,UUS,G,SCAL,DF,A,S,Gl,Ga,B,C,ISC,IC,nSl,Nl,Na, 
. N3,N4,N5) 

IF(IGOTO.EO.O) LOOPCNT— 999 

AliALVSIS ROBULE 

CALL AKALVS 
CDJ"AC3J 

IF {IPLOT.EO.0) GO TO 1031 

L-RITE PLOTTING INFORRATION TO TAPE? 

IF RECUESTED 

J30 IF (ITEH.EQ.SAUITER) GO TO 999 
IF a.KE.l) UaiTE(7,PLOT) 

SAUITcR'ITER 

rag SAvoDJ-osj 

CO 101 n-i,Ni 

£Auxan«x(ii) 

101 CGHTINUE 

CO tea ii-i,Na 
£AUG<II)*G(II > 

102 CONTINUE 

1001 IF (IG0T0.EQ.3) GO TO 1100 
ICC3 CONTINUE 

URITE PLOTTING INFORRATION TO TAPE3 
IF REQUESTED 

1100 IF (IPLOT.EO.0) GO TO 1200 
URITE(7,PL0T) 

NOITER’ITER+1 
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URITEO.IOPT) 

URITEO.ge?) TITLE 
907 FORriAT(4A10) 

URITEO.PLTINPT) 

URITE<6,90€) 

906 FORflATUHl,* PLOTTING INFORNATION URITTEN TO TAPE.*) 
1200 RETURN 
END 



ROUTINE NAME: PLOT 

FILE NAME: PLOTS ( source ) , PLOTB(binary) 

PROGRAM 2 

PROBLEM INDEPENDENT 

OPTION(S): PLOT LINEAR, PLOT CONTOUR, OPTIMIZE WITH PLOTS 

INPUT: TAPES - generated by program one 

Namelist /lOPT/ IPLOT - determines which plotting option 

was selected 

TITLE(1)-TITLE(4) - title for plots 

OUTPUT: OUTPUT (TAPES) - if IPL0T=0, a message "NO PLOTTING" is printed 
COMMON: /PLT/ IPL0T,TITLE(4) 

SUBROUTINE CALL(S): LINPLOT, CONPLOT, HISPLOT 

DESCRIPTION: PLOT is the main routine in program two. It reads in the value 

of IPLOT to determine which plots were requested and executes 
the appropriate subroutine. 
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PTOCIWW PL0T( INPUT, TAPE3. TAPE?. OUTPUT, TAPE5-INPLIT.Ti^(’E 6*0UTPUT) 
COnnON /-PIT/ IPL0T,TITLE(4 1 
MPnELIST /lOPT/ IPLOT 

READ IPLOT TO DETERNIhE UHlCH OPTION 
IS REQUESTED 
READ TITLE FOR PLOTS 

READO.IOPT) 

READ(3,M2) TITLE 
F0RRAT(4A1«) 

CALL APPROPRIATE SUIROUTINE TO PERFORM 
TME PLOTTING 

IF (IPLOT. EQ.l) CALL LINPLOT 
IF (IPL0T.EO.2) CALL CONPLOT 
IF (IPLOT. EQ. 3) CALL HISPLOT 
IF (IPLOT. EQ. 8) URITE(6,9«1) 

Ml FORMAT! 1H1,» NO PLOTTING*) 

STOP 

END 
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ROUTINE NAME: LINPLOT 

FILE NAME: PLOTS( source ) , PLOTB(binary) 

PROGRAM 2 

PROBLEM DEPENDENT: dimension statements 

OPTION(S): PLOT LINEAR 

INPUT: INPUT(TAPE5) - user supplied 

LABELXjLABELY - labels for axes of plot 
TAPE? - generated by program one 
Namelist /PLOT/ NX,XARRAY,YARRAY 
NX - number of data pairs 
XARRAY - the x values of the data pairs 
YARRAY - the y values of the data pairs 

OUTPUT: 0UTPUT(TAPE6) - confirms the creation of the plotting file which 

gets sent to the plotter 

COMMON: /PLT/ IPL0T,TITLE(4) 

SUBROUTINE CALL(S): LRCGOSF Library 

DESCRIPTION: This routine creates the plotting file for the two variable 
linear plots generated by the PLOT LINEAR option. It is 
designed to allow the user to specify any input variable for the 
x-axis and any output variable for the y-axis. 
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SUIROUTINE L INPLOT 

DINENSION XARRAV< ITI1AX4 ). YARRAV t ITNAX4 ) 
COnnON /PLT/ IPL0T,TITLE(4) 

NANELIST /PLOT/ NX.XARRAV, YARRAV 
CALL PSEUDO 

READ DATA 

READ(7,PL0T) 

READ(5,9ei) LABELX, LABELV 
F0RnAT(Aie.5X,Aie> 

NXl-NX+1 

NX2*NX«S 

PLOT 


CALL CALPLT(2.,2..-3) 

CALL ASCALE(XARRAV.8.,NX,l,ie.O) 

CALL ASCALE(VARRAY,8.,NX,1.10.0) 

CALL AXES ( 0 . , d . , 0 • . 8 . , XARRAV ( NX 1 ) , XARRAVt NX2 ) , 1 . 0 , 0 . . LABELX , 
2 -ie,3) 

CALL AXES ( e . , 0 . , 90 . , 8 . . VARR AV ( NX 1 ) « YARRAV ( NX2 ) , 1 . 0, 0 . , LABELV 


2 10,31 

CALL NOTATE < 1 . 75 , 8 . S . . 1 5 , T 1 TLE . 0 . . -40 ) 
CALL L INPLT ( XARRAV, YARRAV , NX , 1 , 1 , 1 , 1 ) 


CALL NFRANE 


.15, 

,.15 



URITE UERIFICATION OF PLOT 


UR1TE(6,905) TITLE 
SMS FORNATdHl.^AlOl 

URITE(6,9eS) LABELX, LABELV 

SMS FORMAT! 1 HO, >TUO UARIABLE PLOT FOR X,A10,t AND >,A10, 
2 tURITTEN TO PLOT FILE*) 

RETURN 

END 
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ROUTINE NAME: CONPLOT 

FILE NAME: PLOTS( source ) , PLOTB(binary) 

PROGRAM 2 

OPTION (S): PLOT CONTOUR 

INPUT: INPUT(TAPE5) - user supplied 

LABELX,LABELY, LABEL! - labels for the x and y axes and definition 

of contours 

TAPE! - generated by program one 

Namelist /PLOT/ XMIN,XMAX,XINC,YMIN,YMAX,YINC,Z,NCONINT 
XMIN - minimum x value 
XMAX - maximum x value 
XINC - increment for x 
YMIN - minimum y value 
YMAX - maximum y value 
YINC - increment for y 

Z - array of values that corresponds to the (x,y) pairs 
NCONINT - the number of contours to plot 

OUTPUT: 0UTPUT(TAPE6) - confirms the creation of the plotting file 

COMMON: /PLT/ IPL0T,TITLE(4) 

SUBROUTINE CALL(S): PL0T3D Library 

DESCRIPTION: Similar to LINPLOT, this routine produces a contour plot of any 

output quantity the user chooses, with respect to any two input 
variables. 
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ooo non 


SUBROUTINE CONPLOT 

REPL Z('«0.40),XY(9),nn(-4) 

COnriON /PLT/ IPL0T.TITLE(4) 

NAHELIST /PLOT/ XMIN.XMAX.XINC, 

CALL PSEUDO 

CALL CALPLT(2.,a..-3) 


VniN.VIIAX.VINC.Z.NCONINT 


READ DATA 


NC«40 

READ(7,PLOT) 

READ(5,903) LABELX, LABELV, LABELZ 
FORnAT(A10,5X,A10,5X,A10) 


PLOT 

(i*<xnAx-xniN)/xiNC+i 
N-(vnAx-vniN)/viNC+i 

L-40 

CALL nAxniN(Z,L.n,N,nn) 

FLO*nn<i)-. 0000001 
Hi-nnca) 

FiNC-cnncei-nna ) )/nconint 
I-FINC 

IF (FINC.GT.1.0.AND.FINC.NE.I) FlNC-I+1.0 
NSET-NHI-0 
ND0T*683 
TSN»0.0 
XVCD'XrilN 
XY(2)*vniM 
XY(3)»Xt1AX 
XY(4)-VHAX 
XY(5)*0. 

XY(6)*0. 

XY(7)*(n-l)/8. 

XY(8)«(N-l)/8. 

XV(9 )* 1 

CALL USC0NTR(Z,L.«.N»FL0.HI,FINC,NSET,NHI,ND0T,XY,TSN, 
IXS.YS) 

c 

C LABEL PLOT UITH PARAMETERS 
C 

IS»3 

LT“3 

UD*1.5 

HD«4. 

AD=-1. 

IC*1 

CALL LABELERCTITLE.NC, IS,LT,UD,HD,RD, AO. AD, IC ) 

UD«1. 

K 0 « 2 . 

AD— 1. 

CALL LABELER! IIHCOHTOURS OF.ll.IS.LT.UD.HD.RD.AO.AD.IC) 
HD-4. 

CALL LAB£LER(LABELZ,NC,IS,L7,UD,HD,RD,A0,AD,1C) 

LT-l 

UO-1. 

HD-4. 

AD--1. 

IC-l 

NC-10 

CALL LABELERCLABELX.NC. IS.LT.UD.HD.RD, AO, AD. IC ) 

LT-4 

UD-1. 

HO-4. 

AO— 1. 

IC-l 

CALL LABELERCLABELY.NC. IS.LT.UD.HD.RD. AO, AD, 1C) 

CALL NFRAHE 

CALL CALPLT(t.,l.,999) 
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ooo 


URITE UERIFICATIOM OF PLOT 


966 

907 


URITE(6.966) TITLE 
F0RFWT(1H1.4A16) 

URITE(6,907) LABELX. LABELY, LABELZ 
F0RI1AT(IH6,»C0NT0UR PL0T»./.4X,»X AXIS • 
. »V AXIS ■ »,A10,/,4X,tCOMTOURS • »,A10) 
RETURN 


END 


»,A10,/,4X, 
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ROUTINE NAME: HISPLOT 

FILE NAME: PLOTS(source) , PLOTB(binary) 

PROGRAM 2 

PROBLEM DEPENDENT ~ dimension statements 

OPTION(S): PLOT CONTOUR 

INPUT: INPUT(^APE5) - user supplied 

NORMOBJ.NORMDV - keys to determine whether or not to normalize 

the objective function and design variable plots. 
TAPES - generated by program one 
Namelist /PLTINPT/ NOITER,ITMAX,NDV,NCON 

NOITER - the number of iterations for the optimization loop 
ITMAX - iteration maximum for CONMIN 
NDV - number of design variables 

NCON - number of constraints 

TAPE? - generated by program one 
Namelist /PLOT/ ITER,OBJ,X,G - this namelist appears NOITER times 

on this file (once for each iteration) 

ITER - iteration number 

OBJ - objective function value 

X - an array of values for the design variables 

G - an array of constraint values 

OUTPUT: OUTPUT (TAPE6) - confirms the creation of the plotting file 

COMMON: /PLT/ IPL0TJITLE(4) 

SUBROUTINE CALL(S): MAXMIN, LINZERO, LRCGOSF Library 
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ROUTINE NAME: 
DESCRIPTION: 


HISPLOT (cont'd) 

This routine plots the history of the optimization process, 
the result from the OPTIMIZATION WITH PLOTS option. The 
plots can be divided into three sets: 

Set 1: Iterations versus Objective Function 

Set 2: Iterations versus Design Variables 

Set 3: Iterations versus Constraints 

Because a design problem may contain both a large number of 
design variables and/or constraints, a maximum of four quantities 
are placed on a single plot in sets 2 and 3. 
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ooo ooo ooo 


SUBROUTINE HISPLOT 

DINENSION SAUX(N1 ).SAUG( N2 ), XARRAVC ITn«X4 ) , VARR«V( ITr!AX4 ) 
DINENSION Xa(ITnAX2,NDU),C2CITrtftX2,hC0N).X2N0Rn(MDU) 
INTEGER SAUITER 
REAL nn<4) 

COnnON /PLT/ 1PL0T,TITLE(4 J 
NAMELIST /PLTINPT/ NOITER, ITf1AX,NDO,NCON 
NAMELIST /PLOT/ SAUITER, SALOBJ, SAUX, SAiJG 
CALL PSEUDO 

READ IN DATA 


READ(S,899) NORHOBJ.NORnDU 
B99 F0RMAT(2I2) 

READ( 3. PLTINPT) 

DO 1000 I-1,N01TER 
RE AD (7, PLOT) 
XARRAV(I)*SAUITER-I 
YARRAV(I)-SAUOBJ 
DO 101 J-l.NDU 
X2(I,J)-SAUX(J) 
lei CONTINUE 

DO 102 J-1,NC0N 
G2(1.J)*SAUG(J} 

102 CONTINUE 
1C08 CONTINUE 

NOITl-NOITER^l 

N0IT2-N0ITER+2 


FIX SCALES FOR X-AXIS 


XARRAYCNOITl )-0.0 

XARRAV(N0IT2)*N0ITER/8. 

XnAJ*8./N0ITER 


OBJECTIUE FUNCTION - NORMALIZE, THEN PLOT 
OBJNORM-YARRAVCi ) 

IF (OBJNORM.EO.0..OR.NORnOBJ.NE.l) GO TO 105 
DO 10S I*1,N0ITER 
YARRAV< I )-YARRAY( I )/0BJN0RM 
leS CONTINUE 
105 CONTINUE 

CALL MAXHINl VARRAY, ITMAX4,N0ITER, 1 ,t1M ) 

CALL ASCALE(nH,8.,a,l,10.0) 

YARRAY(N0ITl)«nn(3) 

VARRAY(N0IT2)-Mn(4) 

CALL CALPLT(2.,2.,-3) 

CALL AXES(0.,0..0..8.,XARRAYCNOIT1 ).XARRAY(NOIT2),XnAJ,0., 

. 10HITERATIONS..15,-10) 

CALL AXES(0.,0.,90.,8-.VARRAY(NOIT1 ),VARRAY(NOIT2),1.0,0., 

. 12H0BJ FUNCTION,. 15, 12) 

CALL LINPLT(XARRAY,VARRAY,N0ITER,1,1. 1,1 ) 

IF (YARRAY(NOITl).NE.O. ) CALL LINZEROCXARRAY.YARRAY, ITMAX4,N0ITER ) 
CALL NOTATE! 1.75. 8. 5,. 15, TITLE. 0., 40) 

IF (OBJNORM.EO.0..OR.NORMOBJ.NE.1) GO TO 107 
CALL NOTATE(9.0,6.75, .15,13HNORnALIZATION.0.,13) 

CALL NOTATE<9.0,6.5,.15,14HFACTOR AND KEY,0.,14) 

CALL NOTATE<9.0,6.0,.15,4HOBJ-.O.,4) 

CALL NUMBER! 9 .65,6.0,.15,OBJNORM,0.,4) 

CALL PNTPLTlll.0,6.08,1,1 ) 

GO TO 108 
107 CONTINUE 

CALL NOTATE(9.0,6.5, .15,3HKEY,0.,3) 

CALL NOTATE(9.0,G,0, .15,3HOBJ,0.,3) 

CALL PNTPLT!9.9,S.08,1,1 ) 

103 CONTINUE 

CALL NFRAME 
URITE(6.900) TITLE 
000 FORnAT!lHl,4A10) 

URITE!6,901) 

901 FORMAT! lH0,tOBJECTIUE FUNCTION*,/,* — *) 

IF !OBJNORM.NE.0..AND.NORMOBJ.EO.l) URITE(e,902) OBJNORM 
S92 FORMAT! IH ,*NORMALIZED BY STARTING UALUE*,E1S.8 ) 

DO 109 K-1,N0ITER,7 
K2-K46 



ooo ooo 


IF (NOITER.LT.KS) K3-N0ITER 
URITE(B,903) ) 

FORrtATdH ,///,» ITER »,8X,I2,6(14X,ia>,/) 

URITE(6,904) 

FORHATdH ,ia0(»-*)/>'.* 95*^^ iL 
URITE(6,905) (YARRAV(KIO,KIC«K,lCa) 

FORtlATdH ,» FUN X.7E16.8) 

CONTINUE 

DESIGN UARIABLES - N0RNALI2E, THEN PLOT 

xaNORnd)-xad,n 
IF (NORNDO.NE.l) GO TO 118 
DO lie j«a.NDu 

IF (Xad,J).EQ.0. ) XaNORNd)«0. 

CONTINUE ^ „ 

IF (xaNORnd).EQ.e. ) go to iie 
DO 111 J-UNDU 
xaNORH(j)-xad.J) 

DO 111 Id.NOITER 
xad.j)«xa(i,J)/x2N0Ri(J) 

CONTINUE 

FIND NAX AND fllN OF ALL ARRAYS TO BE PLOTTED 
CONTINUE 

CALL nAXniN(X2,lTHAX2.N0ITER.NDU.«PI) 

CALL ASCALE(rn,8.,a,1.10-O) 

VARRAY(N0IT1)»HN{3) 

VARRAY(HOITa)-nn(4) 

DO lia I-1,NDU,4 
V-S.S 

DO 113 K-1.4 
I1-I4IC-1 

IF (II.GT.HOV) GO TO 117 
CO 114 J-l.NOITER 
YA!?RAY(J)-Xa(J.in 
CONTINUE 

IF (K.NE.l) GO TO 115 

CALL AXESColfo^fo! .i^.XARRAVCNOITl ),XARRhY(NOITE ) .XFAJ.O. , 

*CALL AXES(0??0T^90y!8??YARRAV(NOITl ) , VARRAY(NOITE ) , 1.0,0., 

. 16HDESIGN VARIABLES, . 15, 16'' 

CALL LIHPLT<XARRAY,VARRAY,hOITER,t,l,K,l ) 

Y-Y-.S 

va-Y+.08 

IF (K.NE.l) GO TO 116 

CALL NOTATEd.75.8.5,.15.TITLE.0.,40) 

IF (XaNORnd).EQ.0..OR.NORnDO.NE.l) GO TO d9 
CALL NOT ATE (9. 0,6. 75, . 15, 13hNORNALIZATION,0. , 13 ) 

CALL NOTATE(9.0:6.5,:15,14HFACTOR AND KEY,©., 14) 

IF^^IxENORNd ).EQ.0. .OR.NORHDO.NE.l ) GO TO 119 
CALL NOTATE(9.0,Y, .15,1HX,0.,1) 

R-II 

CALL NUnBER(9.15,Y, .15,R,0.,-1 ) 

CALL NOTATE(9.3,Y, ,15,1H-,0.,1 ) 

CALL NUNBER(9.5,Y, .15,xaN0RI1< II ),©.,4) 

CALL PNTPLTdl.0,Ya,K,l) 

GO TO 113 

If’*^Ik'!'e0.1 ) CALL NOTATEC9.0,6.5, .1S,3HKEY,0.,3) 

CALL NOTATE(9.0,Y,.15,1HX,0.,1) 

R'll 

CALL NUriBER(9.15,Y,.lS,R,0.,-l) 

CALL PNTPLT(9.7S,Ya,K,l ) 

If”(YARRAY(NO1T1).NE.0. ) CALL LINZER0(XARRAY,VARRAY,ITNAX4,N0ITER) 

CALL NFRANE 

CONTINUE 
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c 

c 

c 


9«6 


SW7 

916 

13S 


99t 


9*9 

122 

121 

C 

c 

c 

c 

c 


214 


215 


213 

217 


212 


C 

C 

C 


911 


912 


913 

222 

221 


URITE RESULTS 
URlTE(6,9d&) 

F0RNAT(1H1,»DESIGN UARIABLES*,/',* ») 

IF (X2HORH(l).EO.0..OR.NORnDU.NE.l) GO TO 135 
URITE (£,987) 

FORHATCIH ,»NORnALlZED BV STARTING UALLIES*) 
URITE(6.9ie) (J,X2NORI1(J),J*l,NDU) 

FORriATClH ,* -I.ElG.B) 

DO 121 K-l.NOITER,? 

K2«K>6 

IF (H0ITER.LT.K2) K2-N0ITER 
URITE(e.983) (KK-1,KK>K.K21 
URITE<6.988) 

FORMATUH ,129(S-<)./,« D UARl) 

DO 122 J«1,NDU 

URITE<6,989 > J, IX21K<,J),KK-K.K2 ) 

FORrMTdH ,S Xt,I2.2X,7E16.8) 

CONTINUE 

CONTINUE 

CONSTRAINTS - PLOT 


FIND HAX AND tllN OF ALL ARRAYS TO BE PLOTTED 

CALL nAXniN(G2, ITHAXa.NOITER.NCON.nri ) 

CALL ASCALE(nn.8.,2,1.10.0) 

VARRAY(N0IT1 )»nn(3) 

VARRAV(N0IT2)-nn(4' 

DO 212 1-1.NC0N.4 
V-6.5 

DO 213 K>1.4 
II-NK-l 

IF (II.GT.NCON) GO TO 217 
DO 214 J-l.NOITER 
VARRAY(J)*G2(J,II ) 

CONTINUE 

IF (K.NE.l) GO TO 215 

CALL AXES ( 0 ! f 0 '. f 0 ! , 8 ! , X ARR AV ( NO I T 1 ) , XARRAV ( NO I T2 ) , XHA J , 0 . , 

, 10HITERATIONS, .15,-10) 

CALL AXES(0..0.,90..8..VARRAV(NOIT1),YARRAY(NOIT2), 1.0,0.. 

. IIHCONSTRAINTS, .15,11 ) 

CALL NOTATE(1.75,8.5, .15,TITLE,0.,40) 

CALL NOTATE(9.0.G.5,.15.3HXEV,0..3) 

CALL L1NPLT(XARRAY,YARRAY,N0ITER,1,1,K,1) 

V-Y-.5 

Y2«Y+.08 

CALL NOTATEt9.0,Y,.15,lHG,0.,l) 

R*ll 

CALL NUNBER(9.15,Y,.15,R,0.,-1) 

CALL PNTPLT(9.75,Y2,K.l ) 

f Nl tP 

IF (YARRAV(NOIT1).NE.0. ) CALL LIN2ER0(XARRAV, VARRAV. ITHAX4.N0ITER ) 

CALL NFRAME 

CONTINUE 

URITE RESULTS 


URITE(6,911) 

FORNAT(lHl,*CONSTRAlNTS»./,» ») 

DO 221 K-1,N0ITER,7 
K2>KT6 

IF (N0ITER.LT.K2) IC2*N0ITER 
URITE(G,903) (KIC-1,KK*K,K2 ) 

URITE<e,912) 

FORHATllH ,120( »-«),/,» CONSTI ) 

DO 222 J>l,NCON 

UR1TE(6,913) J,(G2CKK,J),K<*K.<2) 
FORNATCIH ,t G>,I2,2X,7E16.8) 

CONTINUE 

CONTINUE 

RETURN 

END 
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Additional subroutines accessed by HISPLOT, MAXMIN and LINZERO* are 
listed below. Subroutine MAXMIN finds the maximum and minimum values in a two- 
dimensional array. Subroutine LINZERO draws the zero line on a plot. 


SUBHOUTIME nflxniN(x,HDini,N,NDira,rpt) 
PE«L x(MDinuNDina),nn(<) 
yniN*xnAX*0 
DO 1«91 IM.NDina 
DO 101 J*1,N 

IF (X(J.I ).LT.xniN' xniN-X(J,I) 

IF (X(J,I ).GT.XnAX) XriAX-X(J,I) 

1*1 CONTINUE 
nnci )-xNiN 
nn(a)-xnAx 
RETURN 
END 


SUBROUTINE LINZERO ( XARRAV , YARR AV , flAXa, NO ) 
DIMENSION X ARRAY ( I1AX2 ) , YARRAY ( PIAX2 ) 
DIMENSION X(4),Y(4) 

NOl-NO+1 

NO2-N0+2 

X(l)"Y(l)*Y(a)-0. 

X(a)»XARRAY(NO)+l. 

X(3)-XARRAY(N01 ) 

X(4)-XARRAY(N0a) 

V(3)-YARRAY(N01 ) 

Y(4)»YARRAY(N0a) 

CALL LINE(X.Y,a.l,0,0,0. ) 

RETURN 

END 
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APPENDIX B: 


EXAMPLE 


This appendix contains a sample execution of the OPTIMIZE WITH PLOTS 
option. The example chosen is the Constrained Rosen-Suzuki Function. This 
example is also presented as the first example in the CONMIN manual and was 
chosen for this reason. A listing of all the user defined files is provided, 
followed by the output. 

Analysis module 

c 

c 
c 
c 

c 


SUBROUTINE ANALVS 

CORNON /VARABLE/ X(4) 

COnnON /CONSTRT/ 0BJ.GC3) 

ANALYSIS 

0BJ-X(l)M2-5.»X(l)*X(2)»*3-5.*X{2)*2.»X(3)«a-ai.*X(3)* 

*Gln”(n«2*X(n5x(2)»t2-X(2UX(3)»»2'»X(3)*X(4)»»2-X(4)-8. 

G(2)»X(l)«2-X(l)+a.»X(2)»«2*X(3)»»2*2.*X(4)»»2-X(5)-10. 

G{3)»2.»Xa)»t2*2.»X(l)*X(2)*«2-X«2)+X(3)»*2-X(4)-5. 

RETURN 

END 


Defines module 


SUBROUTINE DEFINE 
URITE(S.l) 

FORNAT( 1H0, 16HDES1GH VARIABLES,/, 

,5H XU)./. 

.5H X(2)./. 

,5H X(3)./, 

.5H X(4)) 

FOR?WT(iHa.l8HOBJECTIVE FUNCTION,/, 

.59H OBJ • XU)«»2 - 5»X(1) ♦ ><^2)«2 

.36H 21»X(3) ♦ X(-4)»*2 ♦ 7IX(4) + 50) 

UR1TE(6,3) 

FORNATdHO.llHCONSTRAINTS,/. 

.SSH GU) - XU)»»2 + XU) ♦ X<2)«2 - 

,3SH ♦ X(3) ♦ X(4)»»2 - XM) - 8,/, 

.59H G(2) • XU)*»2 - XU) ♦ 2»X(2)K2 

,36H ♦ 2»X(4)«»2 - X(4) - 10.<* 

.59H G<3) • 2IXU)»«2 ♦ 2»XU) ♦ , X(2)|»2 - 
.36H - X<4) - 5) 

RETURN 
END 


5<XC2) * 2<X(3)(X2, 


X(2) * 


X(2) ♦ 


X(3)tl!2, 

X(3)XX2, 

X(3)S«2, 
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Input file 


OPTiniZE UlTH PLOTS 
TEST of OPTiniZE UITH PLOTS OPTION 
CCONPAR 

x<a)'i.. 

X<3)*i.« 

Nl-6, 

N2-U, 

ra-«, 

N-t-4, 

rs-8, 

MOU-4, 

HCOM-3, 

KSIDE-#, 

IPRlMT-a, 

NfCG'O. 

riscAL-e, 

DELfUN-0., 

OTAX-ll, 

C£ND 

1 1 >■ End of record 


Job stream 


<D£L1VERV INFOriATION> 


/JOB 

FRANOP.Tieo.cmieooo. 

JSEP,<USER N0.>. 

CHARGE . <CHARGE NO.>.LRC. 

CELryER.<DELIWERY INFOMATION) 

RfLUieOOO) 

REDUCE<-i 

SETHRl-n ^ ^ 

CALUFAOPPRcVSl •EXDEFB, B 1 -EXANAB.Cl *EX1NPT ,UN1 ■ <USER NO . > ) 
REPLACE. PLTIHPT.PLTDATA. 

PLOT.UARIAN 

EXIT. 

/EOR 

SUBANAL(Nl-6,N2*ll) 

PLTLIN(Nl»6,Na» ll, ITHAX'40 ) 

PLTC0N(N1»S,N2»U) 

C?T(N1»6.N2»IX,H3»4,N4»4,N5*8 ) 

/EOR 

HI^L0T(Nl*6!!N2"n»NDU*4,MC0N*3,ITrAX2«4a,ITriAX4«44) 
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FRANOPP OUTPUT 


PREPROCESSOR CONTROL CAROS 


SUBANAHN1-6/N2-11I 

PLUIN(Nl«6»N2»ll»lTHAX«A01 

PLTC0N(N1»6,N2-11> 

0PT(Nl»6»N2»XliN3«4»N4»4»N5«8l 


PREPROCESSOR CONTROL CAROS 
LINPL0T(ITMAX2«42) 

HISPL0T(N1»6/N2»11»NDV-4»NC0N-3»1TMAX«40»ITMAX2»42) 


tn 

to 



cn 


**** 


♦♦♦♦ TEST OF OPTIMIZE WITH PLOTS OPTION 


DESIGN VARIABLES 
X(l) 

X(2) 

X(3) 

X(A) 


OBJECTIVE FUNCTION 

OBJ • X(l)**2 - 5*XU) ♦ X(2)**2 


5*X(2) ♦ 2*XJ3»**2 2X*X(3) ♦ 


X(A)**2 ♦ 7*X(A) ♦ 50 


CONSTRAINTS 

G(l) • X(l)**2 + X(l) ♦ X(2)**2 - X<2) 

G(2) • X(l)**2 - X(l) + 2*X(2)**2 + 

G(3) - 2*X(1)**2 2*X(1) ♦ X(2)**2 - X(2) ♦ 


X(3)*»2 ♦ 

X(3)**2 

X<3)**2 


X(3» ♦ X(A»**2 

♦ 2*X(A)**2 


X(4) - 8 

X(4) - 10 
X(A) - 5 



SCONPAR 


o 

o 
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IE+01> .1E+01» .IE*0l» .1E*01» 0,0, 



'IT 
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* CONNIN * 

* FORTRAN PROGRAN FOR * 

* CONSTRAINED FUNCTION NININIZATION * 


CONSTRAINED FUNCTION MINIMIZATION 
CONTROL PARAMETERS 


IPRINT 

NOV 

2 

4 

LINOBJ 

ITRM 

0 

3 

CT 

-.10000E4-00 


THETA 

.lOOOOE^Ol 

FDCH 

.lOOOOE-01 


ITMAX NCON 
40 3 

N1 N2 

6 11 

CTMIN 

.40000E-02 

PHI 

.SOOOOE^Ol 

FOCHM 

•lOOOOE-01 


NSIOE ICNDIR 
0 5 

N3 N4 

4 4 

CTL 

-•lOOOOE-01 

OELFUN 

.lOOOOE-03 

ALPHAX 

.10000E>00 


NSCAL NFDG 
0 0 

N5 

8 

CTLMIN 

.lOOOOE-02 

OABFUN 

.31000E-01 

ABOBJl 

.lOOOOE^OO 


ALL CONSTRAINTS ARE NON-LINEAR 


INITIAL FUNCTION INFORMATION 

OBJ • ' .310000E>02 

DECISION VARIABLES (X-VECTOR) 

1) .lOOOOE^Ol .lOOOOEtOl .lOOOOEtOl .IOOOOE'^01 


CONSTRAINT VALUES (G-VECTOR) 

1) -.40000E'»01 -.60000E>01 


-.lOOOOE+01 


<ji 



iter • 1 OBJ • .25^B4E+02 

DECISION VARIABLES (X-VECTOR) 

1) «10436E«^01 .10436E+01 .1247964^01 

CONSTRAINT VALUES IG-VECTDR) 

U -.31307E401 -.557B8E401 0. 


ITER ■ 2 OBJ - .12204E402 

DECISION VARIABLES IX-VECTOR) 

1) -.6S498E400 .10325E401 .23572E401 

CONSTRAINT VALUES <G-VECTOR) 

1) 397756*00 -.13275E+01 -.255B0E-12 


iter ■ 3 OBJ ■ .837636*01 

DECISION VARIABLES (X-VECTOR) 

1) .224406*00 .992686*00 . 

CONSTRAINT VALUES IG-VECTOR) 

1) -.113886*01 -.354276*01 -. 


20345E+01 

568436-13 


ITER ■ 4 OBJ > .694206*01 

DECISION VARIABLES (X-VECTOR) 

1) -«34392E*00 .100436*01 .214986*01 

CONSTRAINT VALUES (G-VECTOR) 

1) 0. -.60266E+00 -.21613E-01 


.868476*00 


.138046*00 


-. 318416*00 


-.803886*00 



ITER ■ 5 OBJ • .63271E^01 

DECISION VARIABLES (X-VECTOR) 

1) -.67S66E-01 .10136E«01 .20734E«01 

CONSTRAINT VALUES (G-VECTORI 

1) -.2022SE«00 >.14382E«01 0. 


ITER - 6 OBJ - .61723Et01 

DECISION VARIABLES (X-VECTORI 

I) -.9A581E-01 .99297E^00 .ZOAOOE^OI 

CONSTRAINT VALUES (G-VECTOR) 

1) 0, -.99507E+00 -.53852E-01 


ITER - 7 OBJ ■ .b0706E^01 

DECISION VARIABLES (X-VECTOR) 

1) .74690E-01 .98928E+00 .19478E+01 

constraint values (G-VECTQR) 

IJ -.16766E-01 -.10302E+01 0. 


ITER - 8 OBJ ■ .60218E«01 

DECISION VARIABLES (X-VECTQR) 

1) >.17653E*01 .1003BE«01 .20139E«01 

CONSTRAINT VALUES (6-VECTQRI 

1) -.I7726E-0X -.10338E+01 ,28422E-13 


ITER >9 OBJ ■ .60182E>01 

DECISION variables (X-VECTOR) 

1» .23921E-01 .99428E+00 .19869E+01 

constraint values (G-VECTOR) 

1) -.15891E-01 -.10472E+01 .U747E-02 


-.81323E+00 


-.96346E^00 


-.10562E4^01 


>.97S23E«00 


-.lOlOZEtOl 


a\ 

^ ITER • 10 OBJ • .60133E+01 

DECISION VARIABLES (X-VECTOR) 

1) -.17147E-01 .10055E+01 .20139E+01 - 

CONSTRAINT VALUES IG-VECTORJ 

1) -.15050E-01 -.10270E+01 .29211E-02 


ITER ■ IX OBJ ■ (bOOVSE^Ol 

DECISION variables tX-VECTOR) 

1) .19A41E-01 .99^82E*00 .19908E+01 

constraint VALUES <G-VECTOR) 

1) -.13894E-01 -.1047<(E*01 .34703E-02 


97533E«00 


,10058E+0l 



FINAL QPTIHIZATIOM INFORHATION 


OBJ ■ •600982E4’01 


DECISION VARIABLES IX-VECTDR) 

11 .19A41E-01 .99A82E+00 


.1990IE«01 


CONSTRAINT VALUES IG-VECTOR) 

1) -.13899E-01 -.lOATAE+Ol .3^703E-02 

THERE ARE 2 ACTIVE CONSTRAINTS 
CONSTRAINT NUMBERS ARE 
1 3 


THERE ARE 0 VIOLATED CONSTRAINTS 
TERMINATION CRITERION 

ABS(OBJ(II-OBJII-in LESS THAN DAIFUN FOR 


lOOSBE^Ol 


3 ITERATIONS 


number of ITERATIONS - II 
OBJECTIVE FUNCTION WAS EVALUATED 
CONSTRAINT FUNCTIONS WERE EVALUATED 


78 TINES 
78 tines 



TEST OF OPTIMIZE UITH PLOTS OPTION 
OBJECTIVE FUNCTION 

NORMALIZED BY STARTING VALUE . 31000000E«02 


ITER 

0 

1 

2 

3 

A 5 

6 

OBJ 

FUN 

.lOOOOOOOEtOl 

.eZ109B38E+00 

.391bS39AE*00 

.2bS73788E*00 

.22210A76EA00 .202AA61AE»00 

.197A3AAOEAOO 

ITER 

7 

8 

9 

10 

11 


OBJ 

FUN 

.19A26A30E«00 

.192631B9E+00 

.192SA0ABE«00 

•1923S987E*00 

.193B6»2AE«00 



DESIGN VARIABLES 


NORMALIZED BY STARTING VALUES 
X 1 • .lOOOOOOOE+01 

X 2 • .lOOOOOOOEAOl 

X 3 • .10000000E«01 

X 4 ■ .lOOOOOOOE^OL 


ITER 012 


D VAR 

X 1 .lOOOOOOOE^Ol .10S408S8E401 -.64<i98316E«00 
X 2 .lOOOOOOOEAOl .10436493EA01 .1032469SEA01 
X 3 .lOOOOOOOEtOl .12478812E*01 . 235721T9E AQl 
X 4 .10000000E401 .8b846822EA00 .1380393SE400 


ITER 789 


0 VAR 

X 1 .84640201E-01 -.76S31029E-02 .33920629E-01 
X 2 .98928228E400 .10038076E401 .99427809E AQO 
X 3 .19478042E401 .20138593E«01 . 19868698E AQl 
X 4 -.10S62390E401 -.97522950E400 -.10102266E401 


or 

00 


3 4 9 b 


•23440079E400 

.992679b2E400 

.2034S335E401 

.31840699E400 


10 


.71471336E-02 

•10094562E401 

.20139041E401 

.97933160E400 


-.33391709E+00 
. 1004338 0E401 
•21497900E401 
-.803B8119E400 


11 


.19441181E-01 

.994B1727E400 

•19907692E401 

-.10058298E401 


-.9796980BE-01 

•10136010E401 

.20733989E401 

-.B1323403E400 


-.84981011E-01 

•99246929E400 

.203999B7E*01 

>.96346139E400 



oi 


CONSTRAINTS 


ITER 

0 

X 

2 

3 


9 

6 

CONST 
G 1 
G 2 
G 3 

-.AOOOOOOOE401 

>.60000000E«01 

-.lOOOOOOOEtOl 

-.3098A05AE+01 

-.55673679E+01 

.bA698979E-0I 

-.9007Ab73E«00 

-.13S0Aa06Et01 

-.59993265E-02 

-.U292970E+01 

-.39A81U6E*01 

.29176031E-01 

t32216982E-02 

-.819A3693E«00 

-.1916980AE-0X 

>.19390296E«00 

-.1A994319E«01 

•17A97368E-01 

.6208379BE-02 

-.99686600Et00 

-.37A392A9E>01 

ITER 

7 

8 

9 

10 

11 



CONST 
G 1 
G 2 
G 3 

-.S1735781E-02 

-.1038b957E«01 

.2318S608E-01 

-.79792078E-02 

•.I0SA0228Et01 

.19A9387bE-01 

-.S3X2X6S0E-02 

-.iaS6S98$E«01 

.22331508E-01 

-.92929439E-02 

-.10372328E«01 

.22939169E-0I 

-.13899099E-01 

-.10A73S1AEA01 

•3A703007E-02 




OPriniZE NITH PLOTS OPTION 



ITEAPTIONS 



OPTIMIZE HITH PLOTS OPTION 



ITERATIONS 
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FIGURE 1: SUBROUTINE AND FILE FLOW OF 




ONE 


















OUTPUT 



71 


FIGURE 2: SUBROUTINE AND FILE FLOW OF PROGRAM TOO 
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